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Abstract.  The  1998  Sep  30  and  2000  Apr  4  events  are  both 
large  western  solar  energetic  particle  (SEP)  events,  in  which 
~  1  MeV  proton  intensities  exceed  10^  particles  (cm^  s  sr 
MeV)“^.  Yet  they  exhibit  quite  different  time  variations  of 
elemental  abundance  ratios  and  energy  spectra.  Using  a  model 
of  SEP  transport  through  self-amplified  interplanetary  Alfven 
waves,  we  fit  the  simultaneous  observations  of  protons  and 
ions  of  helium,  oxygen,  and  iron.  We  conclude  that  wave 
amplification  by  high  energy  protons  largely  account  for  the 
differences  in  the  observed  time  histories. 


1  Introduction 

Wave-particle  interaction  has  provided  a  foundation  for  un¬ 
derstanding  the  transport  and  acceleration  of  energetic  charged 
particles  in  astrophysical  plasmas,  from  flaring  loops  on  the 
sun  to  the  interstellar  media.  The  quasilinear  theory  has  long 
been  used  to  deduce  from  interplanetary  (IP)  magnetic  power 
spectra  the  particle  diffusion  coefficients  governing  the  trans¬ 
port  of  solar  energetic  particles  (SEPs)  and  cosmic  rays  in  the 
heliosphere.  Since  the  theory  demands  that  wave  amplifica¬ 
tion  go  hand  in  hand  with  particle  scattering,  we  have  ex¬ 
plored  the  role  of  particle  driven  waves  with  a  simple  model 
of  SEP  transport  coupled  to  IP  Alfven  waves.  We  found 
that  particle  transport  through  self-amplified  waves  accounts 
semi -quantitatively  for  the  complex  yet  remarkably  system¬ 
atic  and  rigidity-dependent  evolution  of  SEP  abundances  in 
the  1998  April  20  event  (Tylka  et  ak,  1999;  Ng  et  ak,  1999a,b). 
Tylka  et  ak  (2000)  suggested  it  also  played  an  important  role 
in  the  evolution  of  particle  spectra  in  that  and  the  1998  Aug 
25  event. 

If  proton-driven  wave  growth  is  important,  its  effect  on 
abundance  variation  and  anisotropies  must  increase  with  pro¬ 
ton  streaming.  This  expectation  is  supported  by  a  compari¬ 
son  of  the  initial  time  dependence  of  SEP  abundances  and 
anisotropy  variations  between  large  and  small  events  (Reames 
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et  ak,  2000,  2001).  More  precisely,  large  streaming  of  > 
10  MeV  protons  is  required  for  wave  amplification  to  affect 
the  transport  of  ~  MeV/amu  ions.  Thus,  the  1998  Sep  30 
and  2000  Apr  4  events  are  especially  interesting  in  that  their 
peak  proton  intensities  are  both  large  at  1  MeV,  but  the  lat¬ 
ter’s  is  smaller  by  a  factor  of  ~  40  at  20  MeV.  In  this  paper, 
we  apply  a  model  of  the  coupled  evolution  of  SEPs  and  IP 
Alfven  waves  to  study  the  different  time  variations  in  the  H, 
He,  O  and  Fe  spectra,  and  H/He  and  Fe/O  ratios  at  a  few 
MeV/amu  in  these  two  events,  with  the  aim  of  assessing  the 
role  of  wave  amplification  and  other  processes. 


2  Description  of  Events 

1998  Sep  30.  The  intensities  of  ~  2.2  MeV  protons  and 
helium  ions,  and  ~  3.7  MeV/amu  oxygen  and  iron  ions  mea¬ 
sured  by  Wind/LEMT  rose  rapidly  to  large  values  (Fig.l, 
top).  The  He/H  ratio  at  ~  2.2  MeV/amu  rose  from  below 
coronal  values  for  6  hrs,  then  decayed.  In  contrast,  Fe/O  at 
~  3.7  MeV/amu  started  well  above  coronal  values  and  apart 
from  a  soft  shoulder,  decayed  for  ~  30  hours  (Fig.l,  bottom). 
In  Fig.  3,  H,  He,  O  and  Fe  spectra  from  500  keV/amu  to  500 
MeV/amu,  compiled  with  data  from  Wind/EPACT/LEMT, 
ACE/SIS,  IMP8/GME,  ACE/EPAM,  ACE/ULEIS  are  shown 
for  three  epochs.  The  associated  W81  flare  occurred  at  1402 
UT  on  Sep  30.  A  shock  arrived  amidst  intensity  decay  at 
~  0700  UT  on  Oct  2  (transit  time  ~  42h). 

2000  Apr  4.  The  2.2  MeV  H  and  He  intensities  rose  more 
slowly  than  in  the  Sep  30  event  to  comparable  magnitudes, 
but  O  and  Fe  intensities  were  much  smaller  (Fig.  2,  top). 
In  contrast  to  the  Sep  30  event,  Fe/O  and  He/H  ratios  both 
decayed  rapidly  from  large  values  and  Fig.  4  showed  much 
softer  particle  spectra  at  3  epochs.  The  associated  W66  flare 
occurred  at  1511  UT  Apr  4,  and  LASCO  observed  a  CME 
with  a  speed  of  ~  1500  km/s  (M.  Andrews,  private  com¬ 
munication).  The  particle  intensities  increased  moderately 
around  shock  crossing  at  1620  UT  Apr  6  (transit  time  also 
~  42h). 
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Fig.  1.  Top:  Symbols  are  Wind/LEMT  observed  H,  He,  O,  and  Fe 
intensities.  Bottom:  Fe/O  and  He/H  ratios.  Curves:  See  Fig.  2 
caption.  1998  Sep  30  1430  UT  is  t=0. 


Fig.  2.  As  in  Fig.l,  for  the  2000  Apr  4  event.  Curves  are  model 
predictions  using  ion  source  with  O/H,  Fe/H,  He/H  of  1,  1,  0.5  x 
coronal  values.  2000  Apr  4  151 1  UT  is  t=0. 


3  Model 


dinates”:  (r,  t)  in  fixed  inertial  frame  and  (yu,  P)  in  the  solar- 
wind  frame, 


SEPs  are  injected  at  a  shock  traveling  through  the  solar  wind 
plasma  with  an  initial  distribution  of  Alfven  waves.  The 
evolution  of  the  phase-space  densities  /s(p,  P,  r,  f)  of  ion 
species  s  and  the  spectral  density  {k,  r,  t)  of  cr-mode  Alfven 
waves  are  coupled  via  quasilinear  theory.  The  model  in¬ 
cludes,  for  particles:  magnetic  focusing,  quasilinear  pitch- 
angle  scattering  by  waves,  solar-wind  convection,  and  adi¬ 
abatic  deceleration;  and  for  waves:  quasilinear  interaction 
with  SEPs,  kinetic  transport  in  the  solar  wind,  and  phenomeno¬ 
logical  damping  of  waves.  Eor  simplicity  we  assume  a  radial 
magnetic  field,  and  constant  solar  wind  and  shock  velocities. 
The  basic  features  of  the  model  have  been  described  in  Ng  et 
al.(1999b). 

The  governing  equations  are  given  below  in  “mixed  coor- 
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Fig.  3.  Spectra  of  H,  He,  O,  and  Fe  in  the  1998  Sep  30  Event  at  3  epochs,  centered  at  17,  33,  and  48  hours  after  onset.  Data  are  compiled  from 
Wind/EPACT/LEMT  (filled  circles),  ACE/SIS  (filled  triangles),  IMP8/GME(filled  squares),  ACE/EPAM  (half-filled  squares),  ACE/ULEIS 
(half-filled  triangles).  The  curves  are  computed  from  the  model. 
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Fig.  4.  As  in  Fig. 3,  for  the  2000  Apr  4  event  at  23,  34,  and  47  hours  after  onset.  IMPS  H  and  He  data  gap  in  panel  1. 


(5) 

^a{k,r,t)  =  2Uk,r,t){kBoV,)/ikoBVA)  (6) 

In  the  above,  t  is  time,  r  heliocentric  distance,  v  ion  veloc¬ 
ity,  P  rigidity,  jj,  pitch-angle  cosine,  (P,  /r,  r,  t)  species- 
s  diffusion  coefficient  in  /r-space,  B  =  BQ{rojr)‘^  magnetic 
held,  k  wavenumber,  p  =  \n{k/B),  Gg  ion  source  term,  7^ 
growth  rate  of  ,  c  light  speed,  e  elementary  charge,  E  pro¬ 
ton  total  energy,  Va  =  VA,oro/r  Alfven  speed,  W  solar-wind 
speed,  ysh  shock  speed,  rsh  =  ro.sh+fVsh  shock  distance  from 
Sun,  =  ±1  for  outward  (inward)  waves,  =  W  +  Va, 
Bo,  ko.  To  and  ro.sh  are  constants,  and  {p,v,  P,  k,  V„)  is 
the  wave-particle  resonance  function  (Ng  &  Reames,  1995). 

In  eq.  (4)  only  /h  appears,  assuming  minor  ions  contribute 
negligibly  to  wave  evolution  (Lee,  1983).  The  parameters 
(^1  <  0  and  (2  <  0  specify  ad  hoc  damping  rates  on  the 
ambient  and  amplihed  waves,  respectively,  and  /?  controls 
the  latter’s  dependence  on  k/B.  The  isotropic  ion  source  G  g 
per  (cm^  MV^  h)  in  eq.  (1)  is  localized  at  the  shock  within 


±Ar/2,  with  Ar  =  0.0125  AU  the  computational  spatial 
grid  size,  such  that 
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The  spectral  index  <5  =  do  +  S'Vsht  is  assumed  identical  for 
all  species.  6  <g;  1  is  the  maximum  seed-particle  fraction  of 
the  solar  wind  Hux  into  the  shock,  n^pro/r^  the  solar-wind 
s-ion  density  at  r  with  ro  =  1  AU,  Pq.s  the  “seed”  ion  rigid¬ 
ity,  rd  the  source  decay  length,  and  E^n  =  Po,kn{ro,sh/rsh) 
the  exponential  rollover  energy  (cf.  the  radial  dependence  of 
maximum  energy  in  Zank  et  al.  (2000)).  These  parameters 
allow  us  to  study  the  effects  of  various  temporal,  spatial,  and 
spectral  characteristics  of  the  shock-accelerated  ion  source. 

For  the  simulation,  we  start  with  empty  ion  distributions 
and  outward  propagating  Alfven  wave  distributions  which 
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are  steady-state  solutions  of  eq.(2),  with  power-law  spec-  harder  than  other  species  at  high  energies.  Tylkaetal.  (2001 
tral  index  a  and  wave  intensity  /qo  at  f=l  AU  and  kjB  =  ICRC)  attribute  similar  features  in  the  2000  July  14  event  to 

2.308  X  10“^  MV  “t  xhe  coupled  equations  are  solved  by  the  shock  acceleration  of  highly  ionized  iron  superthermal 

finite-difference  technique  (Ng  &  Reames,  1994).  left  over  from  previous  impulsive  events. 


4  Comparison  of  theory  with  observation 

We  performed  a  large  number  of  computer  runs  with  the 
model.  For  each  of  the  observed  events,  we  select  the  results 
of  one  run  to  display  as  curves  superposed  on  the  observed 
data  on  the  intensity,  abundance  ratio,  and  spectra  (Figs.  1- 
4).  These  are  not  best  fits,  but  our  preliminary  attempts  to 
compare  theory  and  observation.  We  obtain  reasonably  good 
fit  to  all  observed  histories  of  intensities  and  abundance  ra¬ 
tios  in  Figs.l  and  2,  except  for  the  slow  rise  in  the  observed 
oxygen  intensity  in  the  Apr  4  event.  The  calculated  spectra  fit 
the  observations  on  the  whole,  but  are  softer  than  observed  at 
high  energies,  implying  that  we  could  do  better  with  slightly 
harder  source  spectra.  We  emphasize  that  the  model  predicts 
absolute  intensities  for  all  species;  the  particle  sources  have 
been  specified  with  identical  energy  spectral  shapes,  and  the 
O/H,  Fe/H,  and  He/H  ratios  in  the  source  are  set  to  1,  1,  and 
0.5 X  the  known  coronal  values  (Reames,  1999).  The  proton 
intensity  at  high  energies  critically  determines  the  amplified 
level  of  Alfven  waves  that  scatter  ~  MeV/amu  ions. 

The  following  parameters  are  common  for  both  event  fits: 
nH,o=  5  cm-3,  Po,h=  3.06  MV,  ro,sh=0.1 125  AU,  V;h=0.024  AU  h 
Va,o  =  0.001  AU  h-\  a  =  -1.7,  (2=  -0.30  h“\  and  oxygen 
and  iron  charge  states  of  6.67  and  1 1.67,  respectively. 

For  the  1998  Sep  30  event,  the  other  model  parameters 
are:  h  =  0.003,  (5o  =  -4,  FJo.kn  =  80  MeV,  5'  =  -3  AU“\ 

=  0.2  AU,  Ci=  -0.012  h-\  loo  =  2.5  x  10^  MeV  cm“^ 

P  =  0.25,  and  W  =  0.011  AUh“\  For  the  2000  Apr  4 
events,  these  are:  b  =  0.1,  (5o  =  —5.5,  J7o,kn  =  10  MeV, 

S'  =  -2  A\J-\  u  =  00,  Cl  =  -0.011  h“\  7oo  =  1  X 
10^  MeV  cm“^  /)  =  1,  and  IF  =  0.01  AU  h“\ 

For  the  2000  Apr  4  event,  the  initial  wave  intensities  are  set 
higher  to  fit  the  slower  intensity  rise;  the  source  spectral  in¬ 
dex  (5o  =  —5.5  and  the  e-folding  “knee”  energy  J7o,kn=10  MeV 
result  from  the  observed  soft  spectra.  The  observed  large  pro¬ 
ton  intensity  below  ~1  MeV  then  leads  to  6=0.1,  i.e.  ~  10% 
of  the  solar  wind  flux  is  accelerated  by  the  shock,  if  the 
source  spectra  extend  down  to  5  keV/amu.  Obviously,  the 
source  spectra  must  flatten  below  500  keV/amu,  the  lowest 
energy  in  the  spectra  reported  here.  By  contrast,  b=0.003  for 
the  1998  Sep  30  event  is  quite  reasonable. 

Our  computations  indicate  that  for  the  1998  Sep  30  event 
with  large  intensities  at  >  10  MeV,  unmitigated  proton-driven 
wave  growth  modifies  particle  intensities  more  severely  than 
observed.  Hence  we  introduced  phenomenological  wave  damp¬ 
ing  in  eq.  (4)  to  moderate  wave  growth  and  bring  the  com¬ 
puted  intensities  closer  to  the  observation.  The  fit  parameters 
for  damping  in  the  two  events  are  not  significantly  differ¬ 
ent.  In  both  events,  the  calculated  spectra  are  softer  than 
observed  at  high  energies,  and  the  fits  could  improve  with 
harder  source  spectra.  The  observed  iron  spectra  become 


5  Conclusion 

On  the  whole,  the  model  predictions  agree  well  with  the  ob¬ 
served  evolution  of  the  spectra  and  intensity  of  all  ion  species 
and  the  He/H  and  Fe/O  ratios.  Given  the  wide  range  of  ener¬ 
gies  and  particle  species,  this  overall  agreement  is  encourag¬ 
ing,  especially  since  the  source  abundances  are  chosen  close 
to  the  observed  coronal  composition.  The  success  of  the 
model  in  reproducing  the  contrasting  behavior  of  hard  and 
soft  SEP  events  with  high  ~MeV  proton  intensities  corrob¬ 
orates  the  findings  in  our  previous  studies  contrasting  large 
and  small  SEP  events  (Reames  et  al.,  2000,  2001),  and  lends 
further  support  to  the  role  of  proton  generated  waves. 
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